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Abstract — In this paper we study the controllability of 
networked systems with static network topologies using tools 
from algebraic graph theory. Each agent in the network acts 
in a decentralized fashion by updating its state in accordance 
with a nearest-neighbor averaging rule, known as the consensus 
dynamics. In order to control the system, external control inputs 
are injected into the so called leader nodes, and the influence 
is propagated throughout the network. Our main result is a 
tight topological lower bound on the rank of the controllability 
matrix for such systems with arbitrary network topologies and 
possibly multiple leaders. 

I. Introduction 

Decentralized control of networked multi-agent systems 
has received a considerable amount of attention during the 
last decade. Numerous applications of decentralized con- 
trol laws have been studied including flocking (e.g., [2]), 
alignment and formation conttol (e.g., [l]-[4]), distributed 
estimation (e.g., [6]), sensor coverage (e.g., [5]) and dis- 
tributed control of robotic networks (e.g., [7]), to name a 
few. In a distributed framework, a global task is achieved 
by the local interactions of agents among each other without 
a centralized control. Some tasks only require a common 
agreement between the agents, whereas others may ask for 
agents to achieve some definite configuration in terms of 
their defined states. For such tasks, a fundamental question 
is whether such a decentralized system can be controlled by 
directly manipulating only some of the agents. This question 
motivates our analysis of the controllability of networked 
systems. 

Controllability of networked systems was initially ad- 
dressed in [8], where a connection between the spectral 
properties of the underlying graph modelling a network, 
and the controllability of the system was analyzed. A more 
topological analysis of the problem was later presented in 
[9] with an emphasis on how the symmetry with respect 
to the leader node affects the controllability of the system. 
More general conditions were presented in [10], [11] by in- 
troducing equitable partitions in the analysis. These concepts 
were extended along with additional results in [12]. In [13], 
these equitable partitions were used to obtain an upper bound 
on the rank of the controllability matrix. Recently, distance 
partitions are used in [14] to obtain a lower bound on the 
rank of the controllability matrix for single-leader networks. 

In this paper we analyse leader-follower networks in which 
the agents utilize a nearest-neighbor averaging rule. Some 



agents, called the leaders, support external control inputs 
that ultimately influence the dynamics of all other agents 
namely followers by spreading throughout the network. We 
explore the controllability of the overall system under this 
setting. Our main result is a topological lower bound on the 
rank of the controllability matrix for any graph structure with 
multiple leaders. This lower bound is based on distances of 
nodes to the control nodes, and it can easily be computed 
directly from the network topology, without having to rely on 
any rank test or spectral analysis of the graph. This problem 
was studied for single-leader networks in [14], and a lower 
bound was obtained using the distance partition with respect 
to the leader. In this work, we tackle the general problem with 
possibly multiple leaders by extending the use of distance 
based relationships to such cases. 

The organization of this paper is as follows: Section [XT] 
presents some preliminaries related to the system dynamics 
and algebraic graph theory. In Section [HI] we present our 
controllability analysis. Section [IV] provides an algorithm 
to compute the proposed lower bound on the rank of the 
controllability matrix for arbitrary networks. Finally, Section 
fVl provides the concluding remarks. 

II. Preliminaries 

Consider a networked system of n agents that utilize 
the same nearest neighbor averaging rule, known as the 
consensus equation, to govern their dynamics. For each 
particular agent i, the consensus equation is given as 



(1) 



where xi is the state of agent i, and A/i is the set of 
agents neighboring agent i. Without loss of generality, let us 
assume that Xi € R, and the interactions among the agents 
are encoded via a static undirected graph Q = (V,E). In 
this graph, each node in the node set, V = {1,2, ... ,n}, 
corresponds to a particular agent, and the edge set, E C 
V x V, is the set of unordered pairs (i, j) depicting that the 
nodes i and j are neighbors. In this context, neighbor nodes 
are the ones that have the measurements of each other's 
states. 

The consensus equation provides a simple, yet powerful 
foundation for decentralized control strategies that can be 



utilized in various tasks, including coverage control, con- 
tainment control, distributed filtering, flocking and formation 
control. With all agents utilizing the consensus equation, 
their states asymptotically converge to the stationary mean,if 
and only if the underlying graph is connected [3]. 

Assume that we would like to control this network simply 
by applying external control signals to some of the nodes. 
Without loss of generality, let the first m nodes be the leaders 
taking the external control inputs, and let the remaining (n — 
m) nodes be the followers whose dynamics are governed by 
{l). Let the m dimensional control input be represented by 
vector u. Then, the dynamics of the leader nodes satisfy 



E' 



Xi) + for i = 1, 2 . . . , m. (2) 



where, [u]i denotes the i th entry of the control vector u. 
When the external control signals are applied to the leader 
nodes, their effect on the dynamics propagates to the rest of 
the nodes through the underlying network. 

Our main goal here is to characterize the controllability 
of the overall system under this setting. In particular, we 
are interested in the dimension of the controllable subspace, 
and aim to relate it to the topology of the underlying network 
from a purely graph theoretic perspective. To this end, we use 
some basic tools from algebraic graph theory, in particular 
the degree matrix, the adjacency matrix, and the graph 
Laplacian. 

Let A be the n X n degree matrix associated with the 
graph. The entries of A are given as 



Wi\ if i = j 

otherwise, 



(3) 



where \J\fi\ denotes the cardinality of A/i, and it is equal to 
the number of neighbors of node i. 

The adjacency matrix, A, is an n x n symmetric matrix 
with its entries given as 



l if(i,j)eE 

otherwise. 



(4) 



The graph Laplacian, L, is simply given as the difference 
of the degree and the adjacency matrices, 



L = A - A. 



(5) 



In light of ([T} and (ffj), the dynamics of the leader-follower 
network with m leaders can be given as 



-Lx + Bu, 



(6) 



where x = [x±,X2, ■ ■ ■ , x n ] T is the state vector obtained by 
stacking the states of each individual node, and B is an n x m 
matrix with the following entries 



[B]i 



1 if i = j 
otherwise. 



(7) 



Note that © represents a standard linear time-invariant 
system and it relates the system dynamics to the graph 
topology through the graph Laplacian. 



III. Controllability of Leader-Follower 
Networks 

In this section we will analyse the controllability of the 
system given in ©. In particular, we present relationships 
between the network topology and the rank of the con- 
trollability matrix for such systems. We start this section 
by referring to the results based on the equitable partitions 
presented in [11], [13]. 

A partition of a graph Q = (V, E) is given by a mapping 
7r : V — >■ {C\, C2, . . . , C r }, where denotes the cell that 
node i gets mapped to, and we use dom(n) to denote the 
domain to which tt maps, i.e., dom{n) — {Ci, C2, • ■ • , C r }. 

Definition (External Equitable Partition): A partition tt of a 
graph Q with cells C\, C2, ■ ■ ■ , C r is said to be an external 
equitable partition (EEP) if each node in cell Ci has the 
same number of neighbors in cell Cj for every i ^ j. 

In the controllability analysis, we are particularly inter- 
ested in the maximal leader-invariant EEP of a graph. An 
EEP is said to be leader-invariant if the leader nodes are 
mapped to singleton cells, and such a mapping is said to be 
maximal if no other leader-invariant EEP with fewer number 
of cells exists. Note that for any graph Q there is a unique 
maximal leader-invariant EEP, ir*. Examples of maximal 
leader-invariant EEPs are depicted in Fig. Q] 





Fig. 1. Maximal leader-invariant EEP's for two networks, Q\ and Qi. Q\ 
has a single leader namely 1, whereas Q2 has two leaders namely 1 and 2. 

Maximal leader-invariant EEPs are useful structures in 
the controllability analysis since the states of the nodes that 
appear in the same cell of the maximal leader-invariant EEP 
asymptotically converge to the same value [11]. 

Theorem 3.1 [11] If Q is a connected graph with tt* being 
its maximal leader-invariant EEP, then for all Ci € dom(TT*) 



lim {x k {t) - xi{t)) = 0,Vfc,/ e n* 

t— >oo 



(8) 



In light of Theorem l3.ll one can at most be able to control 
all of the average state values within each cell of the tt*. 
Hence, the cardinality of dom(ir*) provides an upper bound 
on the rank of the controllability matrix as given in [13]. 

Theorem 3.2 [13] Let Q be a connected network, and tt* 
denote its maximal leader-invariant EEP. Given the dynamics 
in (0, the rank of the controllability matrix, T, satisfies 



rank(T) < \tt*\, 
where \ir*\ is the cardinality of dom(7T*). 



(9) 



The upper bound given in Theorem 13.21 is quite useful 
in analyzing the controllability of a leader-follower network. 
For instance, one can conclude that a system is not com- 
pletely controllable if there exists non-singleton cells in its 
maximal leader-invariant EEP. However, all the cells being 
singletons does not necessarily imply that the network is 
completely controllable. 

Next, we present our main result, a lower bound on the 
rank of the controllability matrix when multiple leaders are 
present. In [14], the authors present a lower bound for 
single-leader networks. To this end, they utilize the distance 
partition of an underlying graph with respect to its leader. In 
this partition all the nodes that are at the same distance from 
the leader are mapped into a single cell. It is shown there 
that the rank of the controllability matrix is greater than or 
equal to the number of cells in this partition. 

Theorem 3.3 [14] Let Q be a connected single-leader net- 
work, and ttd denotes its distance partition with respect 
to the leader. Given the dynamics in @, the rank of the 
controllability matrix, T, satisfies 

\n D \ <rank(T), (10) 

where \ttd\ is the number of cells in the distance partition. 

Similar to the single-leader case, the distances of nodes 
from the leaders appear as the fundamental property in our 
analysis. We start our analysis with the following proposition. 

Proposition 3.4 Let Q = (V, E) be a connected network 
with the dynamics in ((6), and let b k be the k th column of 
the input matrix B. Then, for any node i and leader k, 



l(-L) r b k } t = 



if0<r<d ik 
[A r ]ik if r = d ik 



(ID 



where L is the graph Laplacian, A is the adjacency matrix 
of the graph, and di k is the distance of node i to the leader 
node k. 

Proof: 

Using the equality in (0, (—L) r can be expanded as 

r-l 

(-L) r = (A- A) r = A r +J2 (-l) r - m Sm, (12) 

m=0 

where S m is the sum of all matrices that can be represented 
as a multiplication in which A appears m times and A 
appears r — m times. Note that since A and A have only 
non-negative entries, any matrix that can be represented this 
way has only non-negative entries. Moreover, since A is a 
diagonal matrix with positive entries on the main diagonal, 
it doesn't add or remove zeros when multiplied by a matrix. 
Hence, S m has zeros only at the same locations as A rn , and 
the following condition is satisfied: 



Using (Q and ((Till , the i th entry of the vector (A - A) r & fe 
can be expressed as follows: 

[(A-A) r b k ] t = [(A-AY] ik 

r-1 

= [^ r k + ^(-l) r - m [S m k..(14) 

?TI = 

As A is the adjacency matrix of the graph, [-4 r ]ifc is equal to 
the number of paths of length r from node i to node k. Since 
the distance of node i to the leader node k is dik, [A r ]i k — 
for all < r < di k . Hence, ( fT~3b and ([Pil l together imply 
that [{A - AYb k ]i — for all < r < d lk . Furthermore, 
plugging r = di k into (TBI , we get 



[(A - A) dik b k ]i = [A dik ]ik, 



(15) 



where [*4 difc ]ifc is equal to the number of paths with the 
shortest length, dik, from node i to the leader node k, and 
for a connected graph it is non-zero. 

■ 

In a network with m leaders, for each node i we can 
define an m dimensional distance vector, di, that contains 
the distance of node i to each of the leaders as 

iT 



di = [ da d. 



)2 



(16) 



In our controllability analysis, we utilize the sequences of 
these distance vectors, D = (d 1 , d 2 , . . . , d\ D '), where \D\ 
denotes the length of sequence D. In this representation, we 
drop the lower indices corresponding to the node labels, and 
use the super indices to denote the order of the particular 
vector in the sequence. In particular, we are interested in the 
sequences, D, defined by the following rule: 

Rule: For every dP 6 D, there exists an index, k p € 
{1, 2, • • • , m} (where m is the number of leaders), satisfying 



[d% p > [d p ] kp , Vq>p. 
Example: Consider a set of six vectors, 



(17) 
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A vector sequence satisfying the rule in (fTTT i can be 
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For each vector d p in this sequence, the index k p satisfying 
the rule in (TTTb is marked with a circle. Note that here d 1 = 

g and k\ = 1, as the first element of all other vectors 

d q , where q > 1, is greater than the first element of d 1 
which is 0. Similarly, for the second vector in the sequence, 



d 2 = 



we have ki = 2, as the second element of all 



[S„ 







[AT 



0. 



(13) 



the vectors d q for q > 2 are greater than the second element 
of d 2 , and so on. 

Theorem 3.5 Let T> be the set of all distance-vector se- 
quences, D, satisfying the rule given in dTTl ), and \D*\ = 



max \D\ be the maximum length for such sequences. Then 
the rank of the controllability matrix, T, satisfies 

rank{T)>\D*\. (18) 

Proof: 

For a system with n nodes, the controllability matrix is given 
as 

T=[B (-L)B (-L) 2 B ... (—L) n ~ 1 B ] . (19) 

Now, consider vectors of the form 

(-LY*b kp , (20) 

where r p — [d p ] kp , and b kp denotes the k* column of the 
input matrix B. Let d p be the distance vector of node i, 
i.e. d p = di. Then, we have r p = [di]k = di k , and from 
Proposition 13.41 we know that the i th entry of the vector in 
(l20l is non-zero and equal to [A rv ]i kp - Also, for any node j 
with [dj] ki > [di\ki we have the j th entry of this vector equal 
to zero. Using this along with the sequence rule depicted in 
( [T7| i, we conclude that the n x \D*\ matrix 

[(-L)^b kl (-Lpb k2 ... (- L y\°'\b k ^ ], (21) 

has full column rank since each column has a non-zero entry 
that none of the preceding columns have. Note that for every 
p £ {1,2, . . . , \D*\}, we have r p = [d p ] kp < n—l since the 
distance between any two nodes is always smaller than or 
equal to n — l. Hence, each column of the matrix in (EH is 
also a column of T, and rank of T is greater than or equal to 
rank of the matrix in (|2TV Thus, we have rank(T) > \D*\. 

■ 

The lower bound presented in Theorem 13.51 is tight and 
can not be improved for general graphs by only using the 
distances to the leaders. As a rather simple example, let 
us consider a network with a single leader. In that case, 
the distance vectors di are one dimensional, hence the 
longest sequence satisfying the rule in ( fTTT i starts with and 
monotonically increases to the maximum distance from the 
leader. The length of this sequence is equal to the maximum 
distance plus one, which is equal to the number of cells in 
the distance partition with respect to the leader. Thus, for 
one dimensional case this lower bound is equal to the one 
presented in [14]. A couple of examples with multiple leaders 
are depicted in Fig. |2] For those networks, the lower bounds 
on the dimension of the controllable subspaces are computed 
as \D*\ = 5, and = 6, whereas for both systems the 
actual ranks of the controllability matrices are equal to 6. 
Note that in general there is not a unique sequence with the 
maximum possible length, yet we present sample sequences, 
£>* and D|> in FigE] 

By combining the lower bound in Theorem 13.51 and 
the upper bound in Theorem 13.21 we obtain the following 
corollary for the rank of the controllability matrix for any 
connected leader-follower network with the dynamics given 
in ©. 



Corollary 3.6 Let Q = (V, E) be a connected network with 
the dynamics given in ([6]). Let \D*\ be the maximum sequence 
length for distance-vector sequences satisfying the rule in 
( 1771 ), and n* be the maximal leader-invariant EEP of Q. 
Then, the rank of the controllability matrix, T, satisfies 

\D*\ < rank(T) < (22) 

IV. Computing the lower bound 

In this section we present an algorithm to compute the 
lower bound mentioned in the Theorem 13.51 Let S = 
{di, d2, ■ ■ ■ , d n } be the set of all distance-from-leaders vec- 
tors for a given graph. Given these distance vectors, let us 
consider a way of iteratively generating vector sequences 
satisfying the rule in ( TTTb . Let C p be the set of all distance 
vectors that can be assigned as the p th element of such a 
sequence D. According to these definitions, C\ = S. Once 
a vector from C p is assigned as the p th element of the 
sequence, d p , and an index k p satisfying the sequence rule 
is chosen, C p +i can be obtained from C p as 

C p+1 =C p \{deC p \ [d] kp < [d p ] kp }. (23) 

In order to obtain longer sequences, this iteration must be 
continued until C p = 0. However, in general there are too 
many possible sequences that can be obtained this way, and 
it is not feasible to find the maximum sequence length by 
searching among all these possibilities. Instead, we present 
a necessary condition for a sequence satisfying the rule in 
(T% to have the maximum possible length. This necessary 
condition significantly lowers the number of sequences that 
needs to be considered to find the maximum sequence length. 

Proposition 4.1 Let T>* be a maximum length distance 
vector sequence satisfying the rule given in ( 1771 ), then its 
p th entry, d p , satisfies 

[d p ] kp = min [d] kp (24) 

Proof: Assume, for the sake of contradiction, this is 
not true. Then, there exists a distance vector dj £ C p such 
that [dj] k < [d p ] k . Due to the rule of the sequence, dj can 
not be added to this sequence after d p . However, dj can be 
placed right before d p since its index k p satisfies the rule of 
the sequence. Hence, we obtain a longer sequence satisfying 
the rule by placing dj right before d p , which leads to the 
contradiction that D* does not have the maximum possible 
length. ■ 
Note that in obtaining the lower bound, we only care 
about the lengths of sequences, not about their actual entries. 
Hence, if for any di, dj 6 C p we have [di] kp — [dj] kp — 

min fcflfc , then we do not care whether di or di is added to 

dec p p J 

the sequence as d p since the resulting C p +i will be same as 
long 3.S kp is chosen as the index satisfying the rule. Thus, as 
far as the sequence length is concerned, the only important 
decision at each step of the sequence generation is the choice 
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Fig. 2. Leader-follower networks, Q\ and C?2 , each having two leaders namely, 1 and 2. Each node has its 2-dimensional distance vector (distances to the 
leaders) shown next to itself. For both networks, sample maximum length sequences, D\ and D^, satisfying the rule in 1171 are given. For each vector 
in these sequences, the entry corresponding to the index satisfying the rule is circled. 



of k p . Based on this observation, we present an algorithm that 
can be used to compute the lower bound. 

In this algorithm we define a new variable, C, as the set 
of all possible non-empty sets C p that can be obtained at 
step p. Initially this set only includes the set of all the 
distance vectors, S, since there is a unique C\ namely S. 
For each such C p , one can obtain m (number of leaders) 
different C p +i depending on the choice of k p . Once, these 
Cp+i are computed, we remove all the previous C p and store 
the non-empty C p +i sets in C, and continue the iteration. 
Iterations stop when C = 0. We keep a counter variable £ in 
the algorithm and it is incremented by one every time C is 
updated for the next step. Once we reach C = 0, the final 
value of £ gives us the maximum possible sequence length, 
1X3*1. 



Algorithm I 



1 

2 
3 
4 
5 
6 

7 
8 
9 
10 
11 
12 
13 



initialize: C = {S} and £ = 
while C ^ 

C = 

for i = 1 to | C | 
for j — 1 to m 

C(i-i)m+j = Ci \{d £ Ci | [d]j 

end for 
end for 

C = C_\ {C e C I C* = 0} 

c = c 

£ = £ + ! 
end while 
return £ 



min \d] 7 - } 



For instance, consider the network Q\ with two leaders 
shown in Fig. [2] We can represent the flow of Algorithm I 
as a tree structure shown in the Fig. [3] In this tree diagram, 
each node at a given level p corresponds to an element of 



C that is computed in the line 6 of Algorithm I in the p th 
iteration of the while loop. Algorithm will terminate after 
the fifth iteration of the while loop as all those C^s will be 
empty sets. 

V. Conclusion 

In this paper we presented a graph theoretic analysis on 
the controllability of leader-follower networks with possibly 
multiple leaders. In particular, we presented a tight topolog- 
ical lower bound on the rank of the controllability matrix of 
such systems with arbitrary interaction graphs. This lower 
bound is based on the distances of nodes from the leaders. 
We also presented an algorithm to compute this lower bound 
for any leader-follower network. This lower bound may find 
its applications in various problems such as selecting leaders 
in a network that are sufficient to establish a certain level of 
controllability. 
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